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We show that a supersymmetric axion model naturally induces a hybrid inflation with the wa- 
terfall field identified as a Peccei-Quinn scalar. The Peccei-Quinn scale is predicted to be around 
10 15 GeV for reproducing the large-scale density perturbation of the Universe. After the built-in late- 
time entropy-production process, the axion becomes a dark matter candidate. Several cosmological 
implications are discussed. 
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I. INTRODUCTION 

The standard model in particle physics has succeeded 
to describe the physics below the electroweak scale. It is 
not, however, a complete theory because of the theoreti- 
cal problems. One is the strong CP problem. Despite the 
expectation of the existence of the CP-violating #-term in 
the Lagrangian, C = {0g 2 s /32n 2 )G a ^G^ a with 6 ~ 0(1), 
experimentally it is constrained as 6 < 10 -10 . This 
is solved by the Peccei-Quinn (PQ) mechanism [l|, 
in which the axion, pseudo-Nambu-Goldstone boson as- 
sociated with the spontaneous breakdown of the global 
U(1)pq symmetry, dynamically relaxes the 9 to nearly 
zero. Another problem is the gauge hierarchy problem, 
which states that the huge difference between the weak 
scale and the grand unified theory scale requires unnat- 
ural fine tuning. In the framework of supcrsymmetry 
(SUSY) 0], this problem docs not arise. Thus it is well 
motivated that we consider the SUSY axion model. 

On the other hand, cosmological observations revealed 
that the Universe started with inflationary expansion era 
and is now filled with unknown matter, called dark mat- 
ter The cosmological inflation and dark matter can- 
not be accounted for in the framework of the standard 
model, and hence they strongly indicate the physics be- 
yond the standard model. 

In this letter we point out that inflation naturally takes 
place in the SUSY axion model. It takes the form of hy- 
brid inflation where the waterfall field is identified with 
the PQ scalar and the end of inflation is the PQphase 
transition. The idea was already noticed in Ref. @, but 
the discussion there was far from complete, considering 
developments on the SUSY h ybr id inflation model there- 
after [a-flEl]. (See also Refs. [16l [l7j for somewhat simi- 
lar but different models.) Particularly, the PQ symmetry 
breaking scale, /„, is determined as f a ~ 10 15 GeV by the 
condition that the cosmological density perturbation is in 
a correct magnitude. This seems to be too large, since the 
axion coherent oscillation might have abundance much 
larger than the dark matter. However, we point out that 
due to the post inflationary dynamics of the flat direction 
in the scalar potential, the axion is diluted, and it can 
take a role of dominant component of the dark matter. 



Therefore this model provides a simultaneous solution to 
the hierarchy and strong CP problems, inflation and dark 
matter in a simple and unified framework. 



II. MODEL AND COSMOLOGICAL 
IMPLICATIONS 

Our model is described by the following Kahler and 
supcrpotcntial, 



K = |5| 2 + |*| 2 + |*| 2 



(1) 



W = kS{W - fl) + X^XX + kSYY + W , (2) 

where S, W and # arc gauge singlets, X(X) and Y(Y) 
have some gauge charges, and k, X and k are coupling 
constants, which are taken to be real and positive. Here 
we keep minimal Kahler potentials only, and effects of 
non-minimal terms will be discussed later. The constant 
term Wo( = m 3/2-^J> where m 3 / 2 denotes the gravitino 
mass and Mp is the reduced Planck scale) ensures that 
the cosmological constant is nearly zero in the present 
Universe. This supcrpotcntial possesses a global U(1)pq 
symmetry, which is anomalous at the quantum level, and 
also has the U(l)i{ symmetry whose charge assignments 
are shown in Table. HI After $ and $ obtain vacuum 
expectation values (VEV), this PQ symmetry is spon- 
taneously broken and there appears a pseudo-Nambu- 
Goldstone boson, which dynamically cancels the strong 
CP phase and solves the strong CP problem. 

This is nothing other than the SUSY version of the 
hadronic (or KSVZ) axion model [l8|, if X and X have 
color charge. In this case we can choose Y and Y as 
minimal SUSY standard model (MSSM) Higgses : Y = 
H u and Y = Hd- For a certain choice of k, a sizable 
/z-term is generated after S gets a VEV @ , as we will see 
later. 

It is also possible to choose X and X to be MSSM 
Higgses : X = H U ,X = Hd- In this case, the present 
model describes the SUSY version of the DFSZ axion 
model [j~9| . In this case Y and Y are additional chiral 
supcrmultiplets. It is also allowed to introduce some ad- 
ditional chiral supcrmultiplets like heavy quarks in the 
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TABLE I: Charge assignments on the field content. 



KSVZ model. In order to maintain gauge coupling unifi- 
cation, these additional multiplcts may belong to funda- 
mental representations of SU(5). 

One may notice that the first term in the superpo- 
tential introduced to stabilize the PQ scalar at large 
field value coincides with that used for the hybrid in- 
flation [f| 0] , after identification of S with the inflaton 
and ^(^>) with the waterfall fields. Thus we reach the 
interesting possibility : the PQ sector for solving the 
strong CP problem naturally causes inflation. We do not 
need any additional fields and interactions. According 
to a recent analysis including the effect of constant term 
in the superpotential ^ (la ], the correct magnitude of 
the density perturbation is reproduced for f a ~ 10 15 GeV 
and k ~ 10~ 3 if m 3 / 2 — 1 TcV. At first sight this may 
seem to be disappointing, because such large PQ sym- 
metry brea king scale leads to axion overproduction, as is 
well known [2(J[2l[. In this inflationary scenario, the PQ 
symmetry is restored during inflation and broken after 
that. Thus the phase of the axion takes random values 
for different patch of the Universe, and it is not allowed 
to tune the initial misalignment angle to avoid the axion 
overproduction. 

However, the situation is much better than the first 
thought. This is because the late-time entropy produc- 
tion mechanism, which dilutes the axion abundance to 
the acceptable level, is already built in the present model. 
Therefore, the large PQ scale, /„ ~ 10 15 GeV, is rather 
an appealing feature considering that the axion can take 
a role of the dominant component of dark matter after 
the entropy-production process. 

Now we discuss the scalar field dynamics after infla- 
tion. The scalar potential is given by 



V 



I** - fa 
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! |S| 2 (|#| 2 + |#| 2 ). 



(3) 



Here we have taken X = X = since they quickly set- 
tle at the origin due to the Hubble mass term during 
inflation. The global minimum is located at S = and 
= / 2 . In other words, there is a flat direction along 
which the scalar fields do not feel the potential, ensured 
by the U(1)pq symmetry extended to a complex U(l) 
due to the holomorphy. The SUSY breaking effect lifts 
up the flat direction, saxion, and gives a mass of order 
m 3/2 , 



V SB 



CiTO 2 /2 |*| 2 



c 2 m 2 /2 |*| 2 , 



(4) 



where Ci and c 2 are 0(1) constants. This stabilizes the 
flat direction at \^\ ~ ~ f a . We denote deviation 
from this minimum along the flat direction as cr, and call 



it as saxion. The field also receives a finite-temperature 
effective potential, Vp — a 2 T 4 log>I', where a s is the 
QCD gauge coupling constant, coming from two-loop ef- 
fects even if heavy quarks are decoupled from thermal 
bath [22|. 

After inflation ends, the inflaton S and waterfall fields 
^(^>) oscillates around the minimum, S = 0, = |^| = 
f a , noting that the flat direction at this stage obtains a 
mass of k\S\. The scalar degrees perpendicular to this 
direction, which fully mixes with S, decays much ear- 
lier than the saxion since they have masses of ms ~ nf a - 
The decay is induced by the third term in ^ , and the re- 
heating temperature is around Tr ~ 10 11 GcV for ms ~ 
10 12 GcV and k ~ n. After that, the thermal logarithmic 
comes to dominate and drives the saxion to \^>\ ~ a s Mp 
where the effective thermal mass becomes equal to the 
Hubble parameter, and the saxion stops there until the 
thermal effect becomes irrelevant. When the Hubble pa- 
rameter decreases to ~ "13/21 the mass term dominates 
over the thermal correction, and the saxion begins to os- 
cillate around the minimum, \^>\ ~ \^>\ ~ f a , with an 
initial amplitude of cr,; ~ a s Mp. The abundance of the 
saxion coherent oscillation, in terms of the energy density 
to entropy ratio, is then given by 1 
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a s M P 



(5) 



where m a {^ m.3/2) is the saxion mass. This comes to 
dominate the Universe well before the QCD phase tran- 
sition. In the case of KSVZ model, the saxion decays into 
gluons with the rate 2 
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(6) 



Then the Universe is reheated again by the saxion decay. 
The temperature after the saxion decay is estimated as 



T G - 3 McV 



lOTcV 



3/2 



10 15 GcV \ 

fa J' 



(7) 



and hence is compatible with the lower bound on the re- 
heating temperature [26[ for m a > lOTcV. Notice that 
the saxion also decays into a SUSY particle pair and 
then produces the lightest SUSY particles (LSP) nonthcr- 
mally, which easily exceed the dark matter abundance. 



1 Here we have assumed that the S decays before the saxion begins 
to oscillate. Otherwise, the presence of the k\S\ mass term for 
the saxion makes the saxion oscillation amplitude exponentially 
suppressed [23l | . 

2 Decay into two axions, a — > 2a, must be suppressed for successful 
reheating in the KSVZ model. This requires c\ ~ C2 in Eq. l(4j 

MM- 
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Thus we need to introduce small R-parity violation in 
order for the LSP to decay well before BBN begins, or to 
assume SUSY particles are heavy enough not to be pro- 
duced by the saxion decay. In the case of DFSZ model, 
the saxion decays into Higgs pair or fermion pairs. For 
example, the decay width into the lightest Higgs boson 
pair is 



r„ 



J_rr4 

87T ft 



JL 



(8) 



where /i = A(\l/) gives the higgsino mass, and hence we 
obtain 



T a ~ 5 MeV 



m a \ 3 / 2 /10 15 GcV 



ITcV 



fa J \m<r 



(9) 



Thus in this case we need m a ~ 1 TeV and decay into a 
SUSY particle pair can naturally be forbidden. 

Now let us discuss the abundance of the axion, grav- 
itino and axino after the dilution by the saxion decay. 
The axion abundance, in terms of the density parameter, 
after the dilution is estimated as 123 



n n h 2 ~ 5 x 10" 



To \ ( fa 

lMcV J 1 10 15 GeV 



(10) 



hence it is consistent with the WMAP observation of the 
dark matter abundance. This is appealing, since the PQ 
scale f a ~ 10 15 GeV is required for generating the density 
perturbation of the Universe, while this large PQ scale 
leads to the efficient late-time entropy production, mak- 
ing the axion plausible candidate of dark matter. Note 
that the axion does not have an isocurvature perturba- 
tion in this model, since the PQ symmetry is restored 
during inflation. 

As for the gravitino, the y a re produced both ther- 
mally [29j and nonthermally [28l ] from the inflaton decay, 
but diluted sufficiently. The thermally produced grav- 
itino abundance, in terms of the number to entropy ratio, 
is estimated as 



Fs/2 . lxl0 -(15X)' /2 



T R 



T a \ ( a s M P 



10 n GeV I \ lMeV \ a 



(11) 



This satisfies the bound on the unstable gravitino abun- 
dance from BBN m 3/2 Y 3/2 < 10- 13 -10- 9 GeV for m 3/2 ~ 
1-lOTeV [3(1 HJ for an unstable gravitino. 

The axino, which is the fermionic superpartner of 
the axion, might also have significant effects on cosmol- 
ogy I32U331 . The axino abundance from thermal produc- 
tion [34|, after the dilution, is given by 



iq ZlTeV\ 1/2 
Y a ~1 x 10- 19 



/ T a \ / 10 15 GeV \ 

K [lMcVj \ fa J 



10 n GeV 

2 



a H Mr 



a, 



(12) 



In the present model, the axino mass is generated once 
the A-term potential is included : Va = AnSf 2 + h.c. 
with A ~ TO3/2. Then S has a VEV of ~ A/k, and it gives 
an axino mass of m a = k(S) ~ A. Thus the axino mass 
is comparable to the gravitino. If the axino is not the 
LSP, it has a similar lifetime to the saxion in the KSVZ 
model, and it decays before BBN. The constraint is given 
as Y a < 10 -12 so as not to produce too much LSPs. If the 
axino is the LSP, the bound reads m a Y a < 4 x 10~ 10 GeV. 
In both cases, the constraint is satisfied as is seen in 
Eq. pi). 



III. DISCUSSION 

Here we briefly discuss several remaining issues. 

Coleman- Weinberg correction : Since the waterfall 
fields couple to X(X), there arises a Coleman- Weinberg 
(CW) effective potential for the ^ direction [35| . We have 
explicitly checked that this does not modify the (post- 
inflationary scalar dynamics at all. Also the presence 
of Y(Y) affects the potential of the inflaton (S) through 
the CW correction. As long as k is not much larger than 
k, the inflaton dynamics is not much affected. Since we 
need k > k for relaxing to desired vacuum 0. In this 
letter we have assumed k ~ k. 

Stability of heavy quarks : In the hadronic axion model, 
the heavy quarks are stable unless some additional op- 
erators are introduced. Thus the existence of too much 
heavy quarks may be problematic. In our model, how- 
ever, they are not produced in the early Universe effi- 
ciently, since the reheating temperature is lower than the 
heavy quark mass scale. In the DFSZ model, Y(Y) have 
weak scale masses and may also be stable, and are once 
thcrmalizcd after the reheating. But its relic abundance 
is significantly reduced by the saxion decay and no sig- 
nificant cosmological effects arise. 

Non minimal Kdhler potentials : In the hybrid inflation 
model in supergravity, the scalar spectral index n s takes 
a value from 0.98-1.00 for the range of /„ - 10 15 GeV- 
10 16 GeV in the minimal Kahler potential [TU, [l5j]. It is 
possible to make the spectral index more red tilted and 
fall into the best fit range (n s = 0.963 ± 0.012 with 68% 
C.L. Q), by introducing a non-minimal Kahler potential 
K = ks\S\ 4 /Mp and choosing its coefficient as ks ~ 
0.01-0.02. For example, n s = 0.96 (0.95) is obtained for 
f a = 10 15 GeV (4 x 10 15 GeV) for k s = 0.01 [H, [3. 
We can also introduce non-minimal Kahler potentials, 
as K = (/ells'! 2 !*! 2 + /c 2 |S'| 2 |4'| 2 )/M 2 with fa - k 2 - 
O(l). These terms give an additional Hubble mass to the 
saxion. But this does not modify the saxion dynamics at 
all, since the saxion mass is dominated by the k\S\ term 
during the S oscillation, and after S decays these terms 
become irrelevant. 

Domain wall problem : After inflation U(1)pq is bro- 
ken and cosmic strings are formed. After the QCD phase 
transition, domain walls appear which is bounded by the 
strings. Due to the tension of domain walls, strings (and 
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walls) shrink and finally disappear if the color anomaly 
number is one, as in the KSVZ model. This is not the 
case for the DFSZ model, hence we need to introduce 
several heavy quarks in order to make the color anomaly 
number one. 

Axions emitted from strings : The axionic strings con- 
tinuously emit axions with momentum of order of the 
Hubble scale. After the QCD phase transition it obtains 
a mass, and the energy density stored in the axions emit- 
ted by the strings is comparable to that from the coher- 
ent oscillation |36j . Thus these axions may contribute to 
some non-negligible fraction of the dark matter. 

Baryon asymmetry : Because of the inevitable late- 
time entropy production process, any preexisting baryon 
asymmetry is diluted. One possibility to generate the 
baryon asymmetry which survives the huge dilution is to 
rely on the Affleck-Dine (AD) mechanism 37j. Such a 
scenario was studied in detail in Ref. j38j . The baryon- 
to-entropy ratio is expressed as 



n B ^ Scpm 3 / 2 \ipos\ 2 3T g 
s m^cr|^ 4 



(13) 



where Sep ~ C(l) is the effective CP angle, -008 denotes 
the AD field amplitude at the onset of its oscillation, and 
£ is a suppression factor due to the mismatch of the os- 
cillation epochs of the AD field and the saxion. As an 
example, the AD field lifted by the non-renomalizable 
superpotential Wad = ip n /nM n ~ 3 with a cutoff scale 
M corresponds to ?/> os = (H os M n ~ 3 ) 1 ^ n ^ , where H os 
is the Hubble parameter at the onset of AD fields os- 
cillation. In the present parameter, the AD fields be- 



gins to oscillate around T = Tr, and in this case we 
have £ ~ (V"V-Mp/Tr) 3 . Substituting typical values of 
T R - 10 n GeV, (Ji ~ a s M P , M ~ M P and n = 6 yields 
the correct observed value, ub/s ~ 10 -10 . Note that Q- 
balls formed through the AD mechanism decay well be- 
fore BBN and never dominates the Universe because of 
the presence of the saxion. It should also be noted that 
the baryonic isocurvature fluctuation generated by the 
AD mechanism is of the order of ~ (iZ inf /M) 3 / 4 , where 
H- m { ~ nM 2 /Mp is the Hubble scale during inflation, 
and is safely ignored. 

To summarize, we have shown that a simple SUSY 
axion model naturally causes hybrid inflation. The PQ 
symmetry breaking scale is fixed to be around 10 15 GeV 
so that the correct magnitude of cosmological density 
perturbation is reproduced. Taking into account the post 
inflationary dynamics of the scalar fields contained in the 
model, the late-time entropy production necessarily oc- 
curs and the axion coherent oscillation can be the dark 
matter. This provides a solution to the strong CP and 
gauge hierarchy problems and simultaneously explains a 
cosmic inflation and the presence of dark matter. 
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